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Abstract 

The impact of unstable Majorana v T on primordial nucleosynthesis is con- 
sidered. The mass and lifetime of v T are taken in the intervals 0.1-20 MeV 
and 0.001-400 sec respectively. The studied decay modes are v T — > v ^ + <j> and 
v T — > i/ e + 0, where is a massless (or light) scalar. Integro-differential kinetic 
equations are solved numerically without any simplifying assumptions. Our 
results deviate rather strongly from earlier calculations. Depending on mass, 
lifetime, and decay channels of the v T , the number of effective neutrino species 
(found from 4 He), in addition to the 3 standard ones, varies from -2 to +2.5. 
The abundances of 2 H and 7 Li are also calculated. 
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1 Introduction 



Direct experimental data permits the tau-neutrino to be quite heavy, m Vr < 18 
MeV Q, and if the mass indeed lays in the MeV region, the v T would have a 
very strong impact on primordial nucleosynthesis. The recent data from Super 
Kamiokande 0, however, indicates possible oscillations between and either v T 
or a new sterile neutrino, z/ s , with a large mixing angle, sin 26 w 1 and a small mass 
difference 5m 2 = 10~ 2 — eV 2 . If oscillations proceed into v T then both and v T 
are light and the impact of a v T mass on nucleosynthesis is negligible. On the other 
hand it is not excluded that the oscillations go into a sterile v s and, if this is the 
case, v T may be heavy. Of course a heavy tau-neutrino should be unstable because 
otherwise it violates Gerstein-Zeldovich bound 0, and would overdose the universe. 

If the life-time of v T is larger than a few hundred seconds, it can be considered 
as effectively stable on the nucleosynthesis time scale. The role of massive "stable" 
v T in the primordial nucleosynthesis was considered in several papers with chrono- 
logically improving accuracy. In pioneering work M and a subsequent paper M the 
assumptions of kinetic equilibrium for all participating particles and of validity of 
Maxwell-Boltzmann statistics were made (see also the later paper ||), hence the 
problem was reduced to the solution of an ordinary differential equation of Riccatti 
type. However, nonequilibrium corrections to the spectra of both v T and massless 
neutrinos in the case of m Ur in the MeV range happened to be quite significant |7|. |[ 
and a more refined treatment of the problem had to be developed. In ref. |§ the 
simplifying assumption of Maxwell-Boltzmann statistics was dropped in favor of the 
exact Fermi-Dirac one, but it was assumed that kinetic equilibrium is maintained 
for all the species. Nonequilibrium corrections have been treated by one of the au- 
thors pD|| , in the update to B, who found that they do not strongly change their 
original results. Exact numerical solutions of the full system of kinetic equations for 
all neutrino species without any simplifications have been done in refs. [[□], [1^]. In the 
last work a somewhat better precision was achieved and in particular a higher cut-off 
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in particle momenta was used. Also the expressions for some matrix elements of the 
weak interaction reactions with Majorana particles were corrected. It has indeed been 
proven, that nonequilibrium effects are quite significant, almost up to 50%. 

Consideration of the impact of an unstable v T on primordial nucleosynthesis also 
has a long history. In the earlier papers [13-18] some relevant effects were approxi- 
mately estimated, but the accuracy of the calculations was typically rather low. In the 
next generation of papers [19-23] the level of calculations was improved considerably 
but is still subject to some criticism. General shortcomings of these calculations are 
the assumption of kinetic equilibrium for some of the species or the use of Maxwell- 
Boltzmann statistics. The recent paper [23J] makes use of similar simplifications. 
However, the nonequilibrium corrections, as we have seen already for the stable case 
and will see below for unstable v T) are quite significant. The corrections related to 
quantum statistics are typically at the level of 10% [25|. The complete set of kinetic 
equations for the case of the decay v T — > v e + 0, where is a massless scalar was 



numerically solved without any simplifying assumptions in ref. J26| for v T mass in the 
interval 0.1-1 MeV and life-times larger than 0.1 sec. 



In this paper we extend the analysis of ref. |26j to a much larger mass interval, 
0.1 < m UT /MeV < 20, and life-times between 0.001 and 400 sec both for the decays 
of Majorana type tau-neutrino through the channels v T — > v e + and v T — * + 0. 
Our numerical procedure is somewhat more accurate than that of ref. [O], we have 
a larger cut-off in dimensionless momentum y = pa, 20 instead of ~ 13 (a is the 
expansion parameter of the universe), and we correct some expressions for the matrix 
elements squared, as in our previous paper related to the stable v T . As in the preceding 
papers on decaying v T we neglect all other interactions of the scalar 0, except for the 
decay and inverse decay. The exact treatment of the scattering 0z/ <-> 0z/ and the 
annihilation vv <-> 20 presents considerable technical difficulties because of a larger 
number of contributions into the collision integral and because of the non-polynomial 
form of the matrix elements squared of these processes. It will be considered in the 
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subsequent publication |27| . 



The paper is organized as follows. After an overview in section |2] we describe 
in detail our numerical solution in section |3|. The results are presented in section ^ 
followed by a general discussion and conclusion in section || Throughout the paper 
the natural system of units, h = c = k = 1, is used. 

2 Decays of tau neutrino 

We assume that the v T is a Major ana type fermion which is coupled to a scalar 
boson 0, possibly a Majoron ]2SL [2D|], which is light or even massless. The coupling 
of to neutrinos may have diagonal terms as e.g. giv T v T (j) which are important 
for elastic scattering v T + <-> v T + and annihilation v T + v T <-> 20. The non- 
diagonal coupling g a pT v a ( t ) is responsible for the decay of v T into lighter neutrinos, 
v e or (correspondingly a = e or /x). We assume that one of these two couplings 
dominates, i.e. that u T predominantly decays either into v e <j) or v^cf) and consider 
these two possibilities separately. It is implicitly assumed that both v e and are the 
usual active neutrinos. Since chirality is changed by the coupling to a scalar field, the 
corresponding light neutrinos should also be Majorana particles, otherwise new sterile 
states would be produced by the decay. We also assume, that is a weak singlet, 
because the LEP measurements [[3(J of the total decay width of Z° do not leave room 
for any other light weakly interacting particles except for the already known ones. 

There are several possible ways of production of in the primeval plasma. The first 
and evident one is by the decay v T — > + v a . Another possibility is the annihilation 
u T + v T — >• + and the third one is a possible non-thermal production in the course 
of a phase transition similar to the production of axions at the QCD phase transition. 
We neglect the last possibility, assuming that even if (pseudo)goldstone bosons were 
created in the course of the phase transition, the phase transition took place early 
enough so that the created bosons were diluted by a subsequent entropy release in 
the course of the universe cooling down. The rate of 0-production in z/ T -annihilation 
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can be estimated as: 



n,. 



where v is the relative velocity and a ann is the annihilation cross-section. In the limit 
of large energies, s = 4:E^ m 3> rr? Vr it is equal to: a ann ~ (gf/327rs) ln(s/m 2 ) (see 
e.g. ||3l|). One can check that this rate is small in comparison with the universe 



expansion rate H = a /a, if g\ < 10~ 5 . In this case the production of Majorons 
through annihilation can be neglected and they would dominantly be produced by 
the decay of v T . The opposite case of dominant production of 0's by z/ T -annihilation 



and their influence on nucleosynthesis was approximately considered in ref. [|3T . 

The life-time of v T with respect to the decay into massless particles and v a is 
equal to: 

= (2) 

This is faster than the universe expansion rate if T < 4 • 10 10 g a ^/m v , where the tem- 
perature T and m Vr are expressed in MeV. Hence for g a > 10~ 10 the decay is essential, 
i.e. t Vt < H^ 1 , while still T > m Vr . The interval of life-times of u T , that we consider 



below, t Vt = 0.001 - 400 sec, corresponds to g a = (4 • 10" 9 - 6.5 • 10~ 12 ) ym/MeV. 
Thus there is a large range of parameters (coupling constants and masses) for which 
the decay is faster than the expansion while annihilation is effectively frozen. 

3 Numerical solution 

3.1 Kinetic equations 

The basic equations governing the evolution of the distribution functions f a (a = v e , 
u,l, u T , and 0) are discussed in some detail in our previous works The tables 



of matrix elements squared of the relevant weak interaction processes, | A | , are 
presented there. We have some disagreement of our expressions for | A | 2 with those 



of refs. |11|, P6| . In addition to the processes considered for the case of a stable v T 



we should also take decay and inverse decay of v T into account and include an extra 
distribution function, fy. The distributions of photons and electrons/positrons are 
assumed to have the equilibrium form: 

/e i7 = [l±exp(S/T 7 )]- 1 , (3) 

because the electromagnetic interactions are very strong and keep the electromagnetic 
plasma in thermal equilibrium. The behaviour of the unknown function of one variable 
T(t) is found from the law of covariant energy conservation: 

p = -3H(p + P) , (4) 

where the Hubble parameter, H — a/a, is expressed through the total energy density, 
p, as H 2 = (8irp/3m 2 Pl ); here a(t) is the universe expansion factor (scale factor), P is 
the total pressure and m P i = 1.22 ■ 10 19 GeV is the Planck mass. 
The set of integro-differential kinetic equations has the form: 

(dt ~ llpAMPr *) = l T + J t ay > ( 5 ) 
where the collision integral for two-body reactions 1 + 2 — > 3 + 4 is given by the 
expression: 



rscat = _L^- f ^ d3ps d * n 9 1/1 12 

1 OT? 2^, 7-1/0^3 9Z7W0^3 0P./0^13 I 



2E 1 2E 2 (2tt) 3 2E 3 (2ir) 3 2E 4 (2tt) 



12^34 



{2^\ Pl +p 2 - P3 - P 4)F scat (fi, f2, fs, U) , (6) 

F S cat = /3/4(l-/i)(l-/2)-/i/2(l-/3)(l-/4), \A\ 2 is the weak interaction amplitude 
squared summed over spins of all particles, and S is the symmetrization factor which 
includes 1/2 from the averaging over the spin of the first particle (if necessary), 1/2! 
for each pair of identical particles in initial and final states and the factor 2 if there 
are 2 identical particles in the initial state; the summation is done over all possible 
sets of leptons 2, 3, and 4. 

The "decay" part of the collision integrals enters the r.h.s. of the equations for 
fv T -, fu a (a — e or fj), and and has respectively the form: 

j^ecay = / dE^F dec [E VT , E$, E Vt - E4), (7) 

^VrPvr'Vr J \&v T +Pv T )l I 
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jdecay m f°° ^Vv t Vvt 77/771 77. 77. 77 \ 

J u a = Z7 „ ; L a/ , , , - p — *dec{K T , K T - £„ tt > Kj, 



jdecay 



E v a VvT VT J\(my4p„ a )-p Va \ E Ur 

2m r°° dp„ T p Ur 



E^t Vt J\{ 



[, , VV jf VT F dec {E VT , Efr E Vt - E ), (9) 
where m is the mass of z/ r (we omitted the index z/ T ) and: 

F dec {E VTl E h E Va ) = f Vr {E Vr ) [1 + U(E^)} [1 - U a {Eu a )] 

-nmu a {E Va )[\- u t {e Vt )\. (10) 

The contribution of the decay term, I deca v 5 into the collision integral of eq. ([5]) is 
considerably simpler for numerical calculations than the contribution of scattering, 
jscat^ Decause former is only one-dimensional while the scattering terms only can 
be reduced to two-dimensional integrals. 

Four integro-differential kinetic equations for the distribution functions fj(p,t) 
(j = v ti u ei ^) anc ^ e< i- 0) are s °l ve d numerically by the method described in 
refs. |~2], (see also the earlier papers H33|). Technical details of the calculations 
are presented in the following section. 

It is convenient instead of time and momenta to use the following dimensionless 
variables: 

x = ma(t), Dj = Pjd(t), (11) 

where m is an arbitrary parameter with dimension of mass, which we took as m = 
1 MeV, and the scale factor a(t) is normalized so that a(t) = 1/T U = 1/T 7 at high 
temperatures or at early times. In terms of these variables the kinetic equations (||) 
can be rewritten as: 



ii. ,■<),.])(,■.„:) /;"' • /;'" . (12) 
3.2 Initial conditions 

For numerical integration of the kinetic equations we modified our program for stable 
massive tau neutrinos (see the details of the stable case calculation in our paper fl2H ) 
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by including the decay terms eq. (0-11) in equations (|12[) . 

All the calculations of the two-body-reaction part of the collision integral remain 
the same as in our previous work. So we discuss below only the changes to the 
program related to the decay terms. We solve the system of kinetic equations (|T2|) in 
the 'time' interval Xi n < x < x fi n under several assumptions concerning decay. First, 
we assume that the decay was negligible before the initial time Xi n and take the initial 
majoron distribution function equal to zero, f^{xi n ^y) = 0. We can always do this 
for sufficiently small x, because the collision rates in eq. (|T^) are very high due to the 
factor 1/x 2 while the decay terms are proportional to x 2 /t and even for very small 
life-times, r C 1 sec they are suppressed at small x. To determine the proper initial 
moment Xi n we approximately solved the kinetic equation for /^(x,?/) neglecting the 
contribution from (p in the collision integral and assuming equilibrium distributions 
for all other participating particles. Then we can integrate this equation over 
momentum and obtain the following equation for the energy density of majorons at 
early times: 

9 

ydy l + exp( 



d x {x* P<t> )=2C <t> x 2 -^^ — I -f^log — *» ' , (13) 
V^PtotT 2ir 2 Jo ev-l exp(-y) + exp(^^) 



where G$ = 1.221 ■ 6.582^3/(8^) = 2.777 is a constant and x 4 p tot = 10.75vr 2 /30 is 
the total energy density of photons, e , and the three neutrino species in comoving 
frame, and the lifetime, r, is given in seconds. 

For a small value of the product k = (m 2 , x 2 /4) < 1 we can solve this equation 
and obtain: 

P<t> = 2C,£^(^-^)(l-0.08« + 0.006k 2 + 0(k 3 )) 
3x 4 v /p^r 2n 2 24 

w 0.31— (1 -0.08k + 0.006k 2 + 0(k 3 )) . (14) 

XT 



From equation fll4|) we can find the energy density of scalars for a given initial 
time Xi n and compare it with the energy density of equilibrium massless neutrinos 
p eq = 0.57573/x 4 . We will assume that scalars are absent if p^/peq < 0.01. For 



v r lifetime in the range 1 < r/sec < oo we take the initial time x» n = 0.1 for all 
values of the mass. For 0.1 < r/sec < 1 we take Xi n = 0.1 for u T mass in the range 
0.1 < m v jMeV < 2 and x in = 0.05 for 2 < m u jMeV < 20. For r/sec < 0.1 we 
always take X{ n = 0.05. 

3.3 Choice of the final time 

The contributions to the collision integral in the r.h.s. of eq. ( |12"D from two-body 
reactions are suppressed at large x by at least the factor 1/x 2 . We find that at x ~ 50 
all the quantities, which do not feel the decay, reach their asymptotic values so, to be 
on the safe side, we choose the final time for the reactions x sca t = 100. Beyond this 
time all collision terms are taken to be zero. 

The decay continues until all the tau-neutrinos have disappeared. In our runs we 
assume that the decay terminates at Xdecay at which p UT /p eq < 10~ 4 , where p eq = 
0.57573/x 4 is the energy density of one massless equilibrium neutrino. As the final 
time Xfi n we take the maximum of the two, Xdecay and x scat . 

The nucleosynthesis code of ref. |H| requires the final time x ~ 2000 for the 
total energy density, photon temperature and n <-» p rates. We applied a separate 
program, which calculates these quantities, using the energy conservation law and our 
final values for the distribution functions f Vi and the temperature T 7 at time x^ n . 

3.4 Momentum grid 

The particles produced by v T decays have momenta around y = xm UT /2. The decay 
becomes operative and starts to diminish the v T number density at the time t ~ t Ut or 
at xd ~ v^7; if the inverse decay is not important, i.e. if the mass and life-time are 
sufficiently large. Thus for the decay products we should take into account momenta 
up to yo = C^ff^m VT , where the coefficient C depends on the time when the decay 
is completed. 

On the other hand the distribution functions of the tau neutrino and of the mass- 
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less neutrino, which does not participate in the decay, are suppressed at high mo- 
menta. Already at y — 20 we have p(y > 20) / p ~ 10~ 5 for these particles, and we 
therefore choose the momentum cut-off equal to yi = 20. The cases yo > yi and 
yD < yh are considered separately. 

When y D < y L (it corresponds to t Vt < 0.1 sec for m„ T > 1 MeV and t„ t < 1 
sec for m Vr < 1 MeV) we take a logarithmic grid for the Majorons in the region 
0.01 < y < 20 and a linear grid for the neutrinos in the region < y < 20. 

Since the collision integrals for two-body reactions are two-dimensional, their cal- 
culation takes most of the CPU time, and correspondingly we confined ourselves to 
a 100 point grid for the reaction terms. On the other hand, the decay terms are 
one-dimensional and we took a 1600 point grid in their calculation. 

When yn > yi (it corresponds to t Ut > 0.1 sec for m„ T > 1 MeV and t Vt > 1 
sec for vh Vt < 1 MeV) we took a linear grid for the decay products in the interval 
< y < yo while for the neutrino, which does not participate in the decay, the linear 
grid was taken in the range < y < 20. In the interval < y < 20 we took 100 
points in the grid, and expanded with an equally spaced grid with the same dy for 
the decay products up to y = yo- We calculated the reaction terms in the collision 
integrals only in the interval < y < 20, because the decay into modes with y > 20 
becomes important when the contribution of reactions into the collision integral is 
already suppressed by the factor 1/x 4 with x 3> 1 for any mass. 

At y = the analytical expressions for the collision integral for v T differ from those 
at y ^ due to 0/0-uncertainty arising from the factor 1/ E\p\ in front of the integrals 
and similar vanishing factors in the D-functions (see their definition in our previous 
papers [0, |3^|). Because of that, we consider the point y = separately. This permits 
to compare the collision integrals at y = with those in nearby points, y <C 1, in 
order to check that numerical errors are small in the region of small momentum, 
y < 1. This is important for a massive tau-neutrino, because its distribution function 
rapidly changes in the course of evolution. 
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An important check of our program is the comparison of the run with the large 
life-time, t Vt = 400 sec, with that with an infinite life-time, t Vt = oo, (when the 
decay terms in the collision integral are switched off) for masses in the interval 0.1 < 
m^/MeV < 20. The results of both these runs are very close for all masses. 

3.5 Nucleosynthesis 

In order to extract the implications of a massive decaying v T on light element abun- 
dances, we have modified the standard nucleosynthesis code (ref. ||35|| ) in the following 
way. First we import the final photon temperature, which is different from 1.401/x, 
in order to specify the final value of the baryon to photon ratio, 7710 = 10 10 ?t.b/?7, 7 . We 
calculate various quantities at the correct photon temperature and import the values 
to the code. These imported quantities are the total energy density p to t, the 6 weak 
interaction rates for (n p)-reactions, and finally <i(ln a 3 )/c/T 7 (with the account 
of neutrinos) which governs the evolution of the photon temperature. The excessive 
number of equivalent neutrinos, AN, which may be both positive and negative, is 
interpolated from a look-up table as a function of 4 He and 7710- 

4 Results 

The impact on nucleosynthesis of a heavy decaying tau neutrino strongly depends 
upon the decay channel. In the case of the decay into + the most important 
effect is the overall change in the total energy density and the corresponding change 
of the universe cooling rate. Nonequilibrium corrections to the spectra of v e are 
relatively weak for a small life-time, such that practically all v T have already decayed 
at the moment of neutron-proton freezing, T m 0.5 MeV. For a larger life-time some 
nonequilibrium v e would come from annihilation v T + v T — > v e + v e and, as is well 
known, would directly change the frozen n/p-ratio. The distortion of the v e spectrum 
is much stronger in the case of the decay v T — > v e + 0. Correspondingly the results for 
the light element abundances are considerably different for these two different decay 
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channels. 

In figs. 1 the total energy density of all three neutrino species and of the scalar 
is presented as a function of the tau neutrino mass and life-time for asymptotically 
large values of "time" x. It is given in terms of the effective number of massless 
equilibrium neutrinos, N eff = (p UT +p Ur + p Vr + p (/> )/p eq , where p eq = (77r 2 /120)m^/a; 4 
is the energy density of one massless neutrino. To avoid confusion we should stress, 
that this effective number of neutrinos, N e ff, has nothing to do with the excessive 
number of equivalent neutrinos, AN, which is introduced to describe the change in the 
primordial 4 He abundance (see below). The behavior of the "iso- neutrino" curves is 
quite transparent: for a large life-time the relative contribution from nonrelativistic 
u T into p increases as the first power of the scale factor a ~ x. For larger masses 
the v T becomes nonrelativistic earlier, and its contribution into the energy density 
is bigger. On the other hand, very heavy tau-neutrinos annihilate more efficiently 
and thus their frozen number should be smaller than the number density of less 
massive ones. This explains the extrema in fig. la. For small masses and lifetimes 
a plateau is reached, where a fraction of percent uncertainty makes the iso-curves 
wiggle. This uncertainty arises from merging of different programs for different mass- 
lifetime regions. Somewhat less trivial is the fact, that for larger values of mass and 
relatively small life-times the energy density can be smaller than in the standard 
model. This part of the graph is enlarged and presented in fig. lb. A decline in 
the energy density in this range is related to the early decay of u T , when its number 
density is not yet frozen with respect to annihilation and is Boltzmann suppressed. 

In figs. 2 and 3 the time evolution of the energy (a) and number (b) densities 
of different species are presented. The behavior of and <fi reflect the behavior of 
u T . For a small mass and life-time (fig. 2) u T remains relativistic and its number and 
energy densities follow the equilibrium massless ones. For x > 1, when the decay 
starts to operate, some decrease in n Vr is observed due to filling of the previously 
empty scalar states. However, the energy density, p Vr , does not decrease starting from 
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x ~ 5 due to the above mentioned effect of domination of nonrelativistic particles. 
Later, at x > 20, the Boltzmann suppression becomes essential and both n Ur and p Vr 
quickly drop down. For large mass and life-time (fig. 3) the features in the picture 
are more pronounced. The decrease in n Vr and p Vr for 0.1 < x < 0.5 is related to 
the Boltzmann suppression of almost equilibrium v T . After that, the z/ T -annihilation 
becomes frozen and the energy density goes up. Above x = 3 the decay comes into 
play and v T disappear from the primeval plasma. To test the consistency of our 
calculations we checked the conservation of particle number densities and found it 
valid with an accuracy better than 1%. 

Our main results for the decay channel v r — > + are presented in figs. 4. 
We calculated the production of 4 He for different masses and life-times, using the 



nucleosynthesis code modified as described in section |3T5| , and expressed the variation 
of the mass fraction of 4 He in terms of the extra number of neutrino species which 
gives rise to the same amount of 4 He. The results are presented in two graphical 
forms: as a three dimensional graph AiV = F(m VT , t Vt ) (fig. 4a) and as a series 
of curves AiV = F(m UT ) for different values of t Vt (fig. 4b). For large masses and 
low life-times AiV is negative. This is related to the decrease of the energy density 
discussed above. Thus if m Vr = 10 MeV and t Vt = 0.1 sec, the effective number of 
neutrino species at nucleosynthesis is only 2.5. 

Comparison with the results of other groups shows a rather strong deviation. 
We ascribe this to the simplifying approximations made in the earlier papers, which 
might give rise to a significant difference with the exact calculations, and to a better 
accuracy of our numerical calculations, which is typically at the fraction of per cent 
level. For example in the case of m„ T = 14 MeV and r„ T =0.1 sec we obtain for the 
energy density p/p^t = 2.9, while the group |24j] obtained 2.5. In the limit of small 
life-times and masses our results are close to 3.57 (this is the energy density of three 
light neutrinos and one scalar), while the results of are close to 3.9. The effective 
number of neutrino species found from 4 He in our case is 3 + AiV = 2.9 for m Vr = 10 
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MeV and r Vr = 1 sec, while that found in ref. |HJ] is 3.1. The difference is also large 



for m„ T = 10 MeV and t Ut = 0.01 sec: we find 3 + AiV = 2.66 and the authors of pO 
obtained 2.86. 

Consideration of the decay v T — > v e + proceeds essentially along the same lines. 
The energy densities of different particle species in this case are almost the same as in 
the case of the decay v T — > + 0, so we do not present them. However, the influence 
on light element abundances is now very different from the previous case, because 
the nonequilibrium electronic neutrinos produced by the decay strongly shift the n/p 
ratio. In particular an excess (with respect to the equilibrium distribution) of low en- 
ergy electronic neutrinos would create a smaller n/p-ratio and thus a correspondingly 
smaller mass fraction of 4 He. 

The light element abundances for different values of m Vr) t Vt and the baryon-to- 
photon ratio 77x0 = 10 10 ne/n 7 are summarized in figs. 5-11. Massive unstable v T may 
act both ways, diminishing and increasing the mass fraction of 4 He. In terms of the 
effective number of neutrino species it may correspond to a large and positive AN, 
up to AiV = 2.5, or even to the negative AiV = —2 (see fig. 7). As seen on the AiV 
iso-curves (figs. 8), there is a minimum in AiV for lifetime ~ 0.05 sec and mass ~ 4 
MeV. 

Also for this channel we disagree with previous papers. E.g. for m = 0.6 MeV 
and r = 100 sec we find V( 4 He) « 0.244, whereas ref. |§ obtains V( 4 He) « 0.20. 

It was shown in ref. |14| that late decaying v T with t Vt = 10 3 — 10 4 sec and 
ra Vr > 3.6 MeV would distort the deuterium abundance strongly if the decay pro- 
ceed into electronic neutrinos. These u e would create excessive neutrons through the 
reaction v e + p — > n + e + , which would form extra 2 H. This is seen explicitly on 
fig. 9b, where 2 H clearly increases as a function of lifetime. The deuterium produc- 
tion goes quadratically with the baryon density, and we see, that this effect is much 
less pronounced for low 7710- 

In figs. 9-lla the mass fraction of 4 He is presented as a function of m Vr for several 
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values of t Vt for different baryon-to-photon ratio, r] W = 10 10 n£/n 7 = 3,5,7. In 
figs. 9-1 lb, c the similar graphs for the primordial 2 H and 7 Li are presented. 

More graphs showing various elements ( 2 H, 4 He and 7 Li) and as functions of mass 
and lifetime, for both channels u T — > + <fi and v T — > v e + and for 7710 = 1, 3, 5, 7, 9 
can be found on the web-page: |http : / /tac . dk/~sthansen/ decay /| together with 



plots of the n-p reaction rates. 

Finally we changed the initial condition, allowing the majorons to be in equilib- 
rium at Xi n , f^Xin) = l/(exp(y) — 1). This situation corresponds to the majorons 
being produced by some other mechanism prior to the v T decay as discussed in sec- 
tion 2. In this case the inverse decay is more efficient than for frp(xi n ) = and p Vr 
decreases slower as is seen from fig. 12. The change in AN, as compared to the 
case when f^Xin) = 0, varies between 0.4 and 1.0, ANf, =eq = ANf,=o + (0.4 — 1.0), 
depending on mass and lifetime. In particular for long lifetime this difference goes to 
0.57 as expected for all masses. 

5 The role of possible neutrino oscillations 

The considerations of this paper neglect a possible influence of neutrino oscillations 
on primordial nucleosynthesis. If the Super Kamiokande result || is true, then 
is strongly mixed with another neutrino, either v T or a sterile one. If is mixed 
with z/ T , then the mass of the latter is small, m Vr < 160 keV and our results obtained 
for larger masses would be irrelevant. In the other possible case, when oscillations 
are between and a sterile one, u s , the v T mass is not restricted other than by the 
direct laboratory limit |IJ] and our calculations make sense and permit to improve the 
latter significantly. Still the oscillations into a sterile neutrino state would have an 
important impact on nucleosynthesis and would change our results obtained without 
oscillations. We postpone a detailed numerical calculations of the influence of neutrino 
oscillations on primordial nucleosynthesis for the future study and give only some 
simple estimates here (a brief discussion of the impact of the Super Kamiokande 
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results on nucleosynthesis is also given in ref. |34 

In accordance with the measurement |2| the mixing angle is surprisingly large, 
sin 2 26 = 0.8 — 1 and vacuum oscillations are not suppressed at all. Interactions 
with the thermal bath generally result in a suppression of oscillations, except for a 
possible MSW-resonance. This has been studied in a large number of works and a 



list of references can be found in the recent papers |36], [37], |38fl . It is sufficient for our 



purposes to use the simple results of the early papers |40fl , that the oscillations 
v a would be fully developed if sin 4 29\5m 2 \ > 3 • l(T 3 eV 2 |3§ and sin 2 29 ■ 5m 2 > 



t'2 



3-10 _6 eV [Q. Thus for \5m 2 \ < 10~ 3 eV 2 the oscillations z/ M «-> u s , in the case that z/ M 
were produced by the usual weak interaction processes in the primeval plasma, would 
not change our results considerably. The oscillations of muonic neutrinos created in 
the v T decays would not influence nucleosynthesis if the v T life-time is large enough. 
In this limit the oscillations do not change the total energy density of the plasma 
in comparison with the no-oscillation case. The decays with a short life-time and 
the subsequent oscillations would somewhat amplify the production of extra neutrino 
species in the primeval plasma and might be noticeable in nucleosynthesis. 

The oscillations might be very strongly suppressed if the cosmic lepton asymme- 
try (in the case considered it is muon asymmetry) is large in comparison with the 
usually accepted value 10~ 9 — 10~ 10 . The calculations made in this work should be 
correspondingly modified to take into account the change in the distribution func- 
tion, f v , due to the non- vanishing muonic chemical potential. However, even an 
increase of muon asymmetry by 6-7 orders of magnitude would have a negligible ef- 
fect on the kinetics of nucleosynthesis. This is not so for the oscillations of v e «-> u s , 
when a large electric asymmetry might be generated by the oscillations and could 
shift the equilibrium value of n/p-ratio (for a recent work and the list of references 
sec 



36 



If one sterile neutrino species exists, then it is rather natural to expect three sterile 
ones. This possibility could be realized for example if the mass matrix of neutrinos 
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contains both Majorana and Dirac mass terms jy], |4^| . In this case, if v T indeed is 
heavy and unstable, the anomalies in neutrino physics should/could be explained by 
the oscillations in the sector v e — — v s \ — v S 2- Due to these oscillations there might 
be effectively up to two extra neutrino species at nucleosynthesis, but a massive and 
unstable v T which could effectively create a negative AiV would permit to avoid a 
contradiction with the data. 

For the case of neutrino interactions with majorons the contribution into the neu- 
trino refraction index from majorons in the medium would strongly suppress oscilla- 
tions | 43| ) if the z/ T i/ M 0-coupling constant is larger than 10~ 6 — 10~ 7 . It corresponds 
to the very short life-time t Vt ~ 10~ 6 sec (10~ 7 / g) 2 (MeV/m^). 

Thus to summarize, our calculations should be considerably modified if neutrino 
oscillations, as indicated by Super Kamiokande data, indeed exist. However, first, 
it seems premature to make a strong conclusion based on these data and, second, 
there exist regions in the parameter space for which our results with the neglect of 
oscillations remain relevant. 
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Figure Captions: 

Fig. 1 Asymptotic values for the effective number of neutrinos N e ff = (p Ur + 
Pu T + Pu T + P<t>)/ Peq, where p eq = (7n 2 / 120) ttiq/x 4 is the energy density for one massless 
neutrino. In fig. la we present the region 0.1 < r/sec < 400 and 0.1 < m UT /MeV < 20. 
In fig. lb we show the large mass region with decay times in the region 0.01 < r/sec < 
1. In this region the effective number of neutrinos can be less than 3. 

Fig. 2 For m Vr =0.1 MeV and t Ut = 1 sec we present the evolution in time x 
of the energy densities (fig. 2a) and number of particles densities (fig. 2b) in units 
of massless equilibrium neutrinos (p eq = (77r 2 /120)mQ/x 4 and n eq = 0.183mQ/a; 3 ). In 
the time region 1 < x < 6 the massive v T decays to and the majoron 0, but because 
of the small mass, the inverse decay-process delays the complete vanishing of v T until 
time x = 25. 

Fig. 3 For m Vr = 10 MeV and t Ut = 10 sec we present the evolution in time 
x of energy densities (fig. 3a) and number of particles densities (fig. 3b) in units of 
massless equilibrium neutrinos (p eq = (7n 2 / 120) tjIq/x 4 and n eq = 0.183mQ/x 3 ). In 
the time region 0.1 < x < 3 the decay does not play any role, in the region x < 1 
the annihilation of u T into other particles is important. Massive v T decays to and 
majoron <p around x — 3. 

Fig. 4 Relative number of equivalent massless neutrino species, AN = N eqv — 3, 
as a function of v T mass and lifetime r, found from 4 He. In fig. 4a we present a 
2-dimensional surface for AN{m VT ,r). In fig. 4b AN as a function of v T mass for 
different lifetimes including the asymptotic value for stable massive u T (decay time 
r = oo) is presented. 

Fig. 5 4 He for different lifetimes as a function of the v T mass. i] w = 3 (fig. 5a). 
2 H for different lifetimes (fig. 5b). 7 Li for different lifetimes (fig. 5c). 
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Fig. 6 Fig. 6a shows AiV contours found from 4 He as a function of v T mass 
and lifetime. r] 10 = 3. The contours shown correspond to AiV = -0.6, -0.4, -0.2, 
-0.1, 0.0, 0.1, 0.2, 0.4, 0.6. Fig. 6b shows the same as fig. 6a with small lifetimes 
0.001 < r/sec < 0.1. 

u T -> u e + 

Fig. 7a AiV surface found from 4 He as a function of v T mass and lifetime. r] 10 = 3. 

Fig. 7b AiV for different lifetimes as a function of v T mass. 7710 = 3. 

Fig. 8a AiV contours for r] 10 = 3 as a function of v T mass and lifetime. Contours 
follow the increase with the step of 0.2. 

Fig. 8b Same as fig. 8a with small lifetimes 0.001 < r/sec < 0.1. 

Fig. 8c 4 He contours for r] 10 = 3 as a function of v T mass and lifetime. Contours 
follow the increase with the step of 0.05. 

Fig. 9 4 He (9a), 2 H (9b) and 7 Li (9c) as functions of the v T mass for different 
lifetimes. r] W = 3. 

Fig. 10 4 He (10a), 2 H (10b) and 7 Li (10c) as functions of the v T mass for different 
lifetimes. 7710 = 5. 

Fig. 11 4 He (11a), 2 H (lib) and 7 Li (11c) as functions of the v T mass for different 
lifetimes. 7710 = 7. 

Fig. 12 Energy density of u T and <fi as functions of time x for the cases with 
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Figure la. 
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